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Abstract. We propose an analytical method to study the entangled spatial and spin dynamics of interacting 
bimodal Bose-Einstein condensates. We show that at particular times during the evolution spatial and spin 
dynamics disentangle and the spin squeezing can be predicted by a simple two-mode model. We calculate 
the maximum spin squeezing achievable in experimentally relevant situations with Sodium or Rubidium 
bimodal condensates, including the effect of the dynamics and of one, two and three-body losses. 

PACS. PACS-03.75.Gg Entanglement and decoherence in Bose-Einstein condensates - PACS-42.50.Dv 
Quantum state engineering and measurements - PACS-03.75.Kk Dynamic properties of condensates; col- 
lective and hydrodynamic excitations, superfluid flow - PACS-03.75.Mn Multicomponent condensates; 
spinor condensates 



1 Introduction 

In atomic systems effective spins are collective variables 
that can be defined in terms of orthogonal bosonic modes. 
In this paper the two modes we consider are two different 
internal states of the atoms in a bimodal Bose-Einstein 
,5^ ' condensate. States with a large first order coherence be- 
^^ , tween the two modes, that is with a large mean value 
^ • of the effective spin component in the equatorial plane 
Ci I of the Bloch sphere, can still differ by their spin fluctua- 
tions. For an uncorrelated ensemble of atoms, the quan- 
tum noise is evenly distributed among the spin compo- 
nents orthogonal to the mean spin. However quantum cor- 
relations can redistribute this noise and reduce the vari- 
ance of one spin quadrature with respect to the uncorre- 
lated case, achieving spin squeezing [HE]- Spin-squeezed 
states are multi-particle entangled states that have prac- 
tical interest in atom interferometry, and high precision 
spectroscopy [3j. Quantum entanglement to improve the 
precision of spectroscopic measurements has already been 
used with trapped ions [4] and it could be used in atomic 
clocks where the standard quantum limit has already been 
reached [5]. 

A promising all-atomic route to create spin squeez- 
ing in bimodal condensates, proposed in [5], relies on the 
Kerr-type non linearity due to elastic interactions between 
atoms. Quite analogously to what happens to a coherent 
state in a nonlinear Kerr medium in optics [7], an initial 
"phase state" or coherent spin state, where all the effective 



spins point at the same direction, dynamically evolves into 
a correlated spin-squeezed state. A straightforward way to 
produce the initial phase state in a bimodal condensate is 
to start with one atomic condensate in a given internal 
state a and perform a 7r/2-pulse coupling coherently the 
internal state a to a second internal state b [8j. However, as 
the strength of the interactions between two atoms a — a, 
b — b and a — b are in general different, the change in the 
mean field energy excites the spatial dynamics of the con- 
densate wave functions. In the evolution subsequent to the 
pulse, the spin dynamics creating squeezing and the spa- 
tial dynamics are entangled [6,9,, 10,11 and occur on the 
same time scale set by an effective interaction parameter 
X- This makes it a priori more difficult to obtain simple 
analytical results. 

In this paper we develop a simple formalism which al- 
lows us to calculate analytically or semi analytically the 
effect of the spatial dynamics on spin squeezing. In Sec- 
tion[2]we present our dynamic model. Using our treatment 
we show that at particular times in the evolution the spa- 
tial dynamics and the spin dynamics disentangle and the 
dynamical model gives the same results as a simple two- 
mode model. We also identify configurations of parameters 
in which the simple two mode-model is a good approxi- 
mation at all times. Restricting to a two-mode model, in 
Section[3]we generalize our analytical results of [I2j on op- 
timal spin squeezing in presence of particle losses to the 
case of overlapping and non-symmetric condensates. 
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In Sections m and [5l we apply our treatment to cases of 
practical interest. We first consider a bimodal ^^Rb con- 
densate. Rb is one of the most common atoms in BEC 
experiments and it is a good candidate for atomic clocks 
using trapped atoms on a chip [T3]. Restricting to states 
which are equally affected by a magnetic field to first or- 
der, the most common choices are \F = l,m = —I) and 
|i^ = 2,771 = I) which can be magnetically trapped, or 
\F = l,m = 1) and \F ~ 2,171 ^ —I) that must be trapped 
optically but for which there exists a low-field Feshbach 
resonance which can be used to reduce the inter-species 
scattering length [l4l[T5]. Indeed a particular feature of 
these Rb states is that the three s-wave scattering lengths 
characterizing interactions between a — a, b — b and a — b 
atoms are very close to each other. A consequence is that 
the squeezing dynamics is very slow when the two conden- 
sates overlap. The inter-species Feshbach resonance can be 
used to overcome this problem and speed up the dynamics 

In schemes involving the \F = 2,m = ±1) of Rubid- 
ium, the main limit to the maximum squeezing achievable 
is set by the large two-body losses rate in these states. As 
a second case of experimental interest we then consider 
Na atoms in the \F = l,mp = ±1) states |6j. Although 
theses states have opposite shifts in a magnetic field, they 
present the advantage of negligible two-body losses. Using 
our analytical optimization procedure, we calculate the 
maximum squeezing achievable in this system including 
the effect of spatial dynamics and particle losses. 

In Section [5] we examine a different scenario for Rb 
condensates in which, instead of changing the scattering 
length, one would spatially separate the two condensates 
after the mixing 7r/2 pulse and hold them separately dur- 
ing a well chosen squeezing time. An interesting feature 
of this scheme is that the squeezing dynamics acts only 
when the clouds are spatially separated and it freezes out 
when the two clouds are put back together so that one 
could prepare a spin squeezed state and then keep it for 
a certain time [12]. State-selective potentials for *^Rb in 
\F = l,m = -1) and |F = 2,m = 1) [H] have recently 
been implemented on an atom chip, and such scheme could 
be of experimental interest. 



2 Dynamical spin squeezing model 

In this section we develop and compare dynamical models 
for spin squeezing. No losses will be taken into account in 
this section. 



where h^ is the one-body hamiltonian including kinetic 
energy and external trapping potential 
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The interactions constants g^e' are related to the corre- 
sponding s-wave scattering lengths g^e' = 'iirh^aee' /M 
characterizing a cold collision between an atom in state 
e with an atom in state e' {e,e' = a,b), and M is the 
mass of one atom. 

We assume that we start from a condensate with N 
atoms in the internal state a; the stationary wave function 
of the condensate is (j)o{r). After a tt/2 pulse, a phase state 
is created, which is our initial state: 
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where Cq, Ct are mixing coefficients with |Cap-|-|Cbp — 1 
and the operator al, > creates a particle in the internal 
state a with wave function 0o- To describe the entangled 
evolution of the spin dynamics and the external dynamics 
of the wave functions, it is convenient to introduce Fock 
states with well defined number of particles in \a) and 
1 5), these numbers being preserved during time evolution 
subsequent to the mixing pulse. Expanded over the Fock 
states, the initial state ^ reads: 
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where Nb — N ~ Na, and 
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(5) 
Within an Hartee-Fock type ansatz for the A^-body 
state vector, we calculate the evolution of each Fock state 
in gl). We get [9]: 



\Na ■.(t)o,Nb:(t)a) 



^-iA{Na,Ny,t)/h 



x\Na : MNa,Nb;t),Nb : MNa,Nb;t)) , (6) 

where (f>a{Na,Nb;t) and 4>b{Na,Nb;t) are solutions of the 
coupled Gross-Pitaevskii equations: 



ihdtCl), = [h, + {N, ~ l)geMe 

here with the initial conditions 



Kg,A't>, 



Pl'/'e (7) 



2.1 State evolution 

We consider the model Hamiltonian 

e—a,b 
+ 9ab-ikWJja1pb (1) 



</'a(0) = ^b{0) = 00 , (8) 

and the time dependent phase factor A solves: 



d 



A{Na,Nb;t) = - E ^e{N, - 1)^ I d^Acj,, 



NaNbgab [ d^r\cba\''\(l)b\'' ■ (9) 
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With this treatment we fuUy include the quantum dynam- 
ics of the two condensate modes a and 6, as one does for 
the simple two modes model, but also including the spa- 
tial dynamics of the two modes and their dependence on 
the number of particles. The approximation we make is 
to neglect all the other modes orthogonal to the conden- 
sates which would be populated thermally. An alternative 
method is to use a number conserving Bogoliubov theory 
that explicitly includes the operators of the condensates 
as in p^. In that case all the modes are present but the 
modes orthogonal to the condensates are treated in a lin- 
earized way. In [lOj , the author compares the number con- 
serving Bogoliubov approach to our approach using many 
Gross-Piaevskii equations, also used in [6], and he finds 
very similar result for the spin squeezing. He also finds 
that within the Bogoliuobov approximation the thermally 
excited modes strictly do not affect the squeezing in the 
scheme we consider here. If the number conserving Bogoli- 
ubov has the advantage of being systematic, our approach, 
supplemented with a further approximation (the modulus- 
phase approximation introduced in Sect. 12. 3| ) allows us to 
get some insight and obtain simple analytical results. 




x^ 




2.2 Calculation of spin squeezing 

The effective collective spin of a two-components BEG can 
be represented on the Bloch sphere as shown in FigdjTop). 
Formally, we introduce three spin operators in terms of 
field operators [6] 



S. 



S, 



-Jd^r[^jl{r)i:a{r)+i^i{r)Mr)], 
- J d^r[i;l{r )ijja{r ) - i^Kr )Mr )], 
- / d3r[V;t (r )i;air ) - 4'lir)Mr)]- 



(10) 

(11) 
(12) 



Definitions p0|) - p2)l explicitly take into account the spa- 
tial wave functions of the condensate and depend in par- 
ticular on the overlap between the two modes. 

Referring to the FigmTop) we introduce the polar an- 
gles d and ip giving the direction z' of the mean spin; i? de- 
termines the relative mean atom number in the two inter- 
nal states, cos^ — |Cap — IC^p, while the azimuthal angle 
ip corresponds to the relative phase between the compo- 
nents, ip = arg(C*Cf,). 

The minimal variance of the spin in the plane (x',y') 
orthogonal to the mean spin AS'^ ^j^^, represented in Fig. 
[TjBottom), is given by 



ASl 
where we introduced 



,rnin n ^ x' 



ASl,-\AS^^\) 



D— — Jx' 



iSy' . 



(13) 



(14) 



The degree of squeezing is then quantified by the pa- 
rameter 6,3] 



Fig. 1. (Top) average spin. (Bottom) variance of the the spin 
components in the plane orthogonal to the mean spin. 



where (S) is the length of the average spin. 

When expressed in the original frame of reference, the 
minimal variance in the orthogonal plane is: 

^^Lmin = ^(cos' ^ cos^ p + sin^ ^)AS^ 

+ -r(cos^ d sin^ ip + cos^ 'p)AS?, 
2 " 

+ - sin^ i^AS"^ sin^ ■& sin 2ipAxy 



- sin 2^9 cos ipA^ 
-VA'^ + B^ 



— sin 2'd sin ip Ay z 



(16) 



where 

A = (sin^ ip - cos^ -d cos^ v)ASl 

+ (cos^ ip - cos^ d sin^ ip)ASl " ^i^^ '^^'^z 

(1-1- cos^ z9) sin 2(pAxy H — sin 2-3 cos (pA^ 

1 



— sin 2i? sin ipAy^ ; 



(17) 



B = cost} sin 2p{ASx ~ ASy) - cos ■& cos 2ipAxy 
— sin "(9 sin ipAzx + sin i!) cos ipAyz ; (18) 

and where we introduced the correlations 

A, - {S,S, + S,S,) - 2{S,) (S,), i^j=x,y,z. (19) 

The spin squeezing is then calculated in terms of averages 
of field operators products, with the state of the system 
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at time t, obtained by evolving equation (j4]) with equation 
([6]). To calculate the averages one needs to compute the 
action of the field operators ipa i^b on the Fock states ^ 



Mr) \Na ■■ MNa,m),Nb : MNa,^)) 

= MNa,Nb,r)^a 

X \Na~l:MNa,Nb),Nb:MNa,Nb)), (20) 

Mr) \Na ■■ MNa,Nt,),Nb : MNa,Nb)) 

= MNa,Ni„r)y/N'b 

X \N,:MNa,Nb),Nb-l:MNa,Nb)). (21) 

The explicit expressions of the averages needed to calcu- 
late the spin squeezing parameter are given in Appendix 
[K\ These quantum averages correspond to an initial state 
with a well-defined number of particles N. In case of fluc- 
tuations in the total number of particles where the density 
matrix of the system is a statistical mixture of states with 
a different number of particles, a further averaging of N 
over a probability distribution P{N) is needed pilTT]. 

2.3 Dynamical modulus-phase approach 

In principle, equations (O-® can be solved numerically 
for each Fock state in the sum equation (j4|) , and the squeez- 
ing can be computed as explained in the previous section. 
However, for a large number of atoms and especially in 
three dimensions and in the absence of particular symme- 
tries (e.g. spherical symmetry) this can be a very heavy 
numerical task. To overcome this difficulty, in order to 
develop an analytical approach, we can exploit the fact 
that for large N in the initial state ^ the distributions of 
the number of atoms iV^ and Ni, are very peaked around 
their average values with a typical width of order VN. 
Moreover, assuming that possible fluctuations in the total 
number of particles are described by a distribution P{N) 
having a width much smaller than the average of the total 
number of particles TV, we can limit to Na and Nb close 
to Na = \Ca\^N and Nb = IM^N- We then split the 
condensate wave function into modulus and phase 



I exp(i6'e) 



e — a,b , 



(22) 



and we assume that the variation of the modulus over the 
distribution of N^ can be neglected while we approximate 
the variation of the phase by a linear expansion around 
Ng [9] . The approximate condensate wave functions read 



(j>e{Na,Nb) ~ ^eCXp 



J2 (M ^ NMdNjs)N.,N, 



s'—a,b 



(23) 
where (j), = M^a ^Na,Nb = Nb). 

The modulus phase approximation takes into account, 
in an approximate way, the dependence of the condensate 
wave functions on the number of particles. It is precisely 



this effect that is responsible of entanglement between spa- 
tial dynamics and spin dynamics. 

As explained in Appendix |B1 all the relevant averages 
needed to calculate spin squeezing can then be expressed 
in terms of (p^ and of three time and position dependent 
quantities: 



Xdir) - - [{dN^ - dNjiOa - 0b)M,N^ , (24) 

Xs{r) = \ [{dN^ + dN,){ea - efc)]^„,^. - (25) 
Xo(r) = i [{d^^ - dNjiea + Ob)]j^^^^^ . (26) 



In some cases (see Sect. 12. 4p these quantities can be ex- 
plicitly calculated analytically. To calculate the squeezing 
in the general case, it is sufficient to evolve a few coupled 
Gross-Pitaevskii equations ([7]) for different values of Na, 
Nb, to calculate numerically the derivatives of the phases 
appearing in ([M)l -(|26 p . Although we do not expect a per- 
fect quantitative agreement with the full numerical model 
for all values of parameters, we will see that the analytical 
model catches the main features and allows us to interpret 
simply the results. 

In the particular case of stationary wave functions of 
the condensates, the parameters Xd, Xs and xo become 
space-independent : 



St _ [(^Jv„ - dN,)M - M]n^,N, . 
Xd — :^ ^ 



2ft 

[{On^ +dNjM -Mb)] 



Na,N^ 



2h 



xf = xf . 



(27) 

(28) 
(29) 



In this case we recover a simple two-mode model. Equa- 
tions ((77|) - (P5)) will be used in section [31 In that contest 
we will rename x^*/i = — x and x|*/t = — x to shorten the 
notations. 

To test our modulus-phase dynamical model, in Fig. [21 
we consider a situation in which the external dynamics is 
significantly excited after the 7r/2 pulse which populates 
the state b. Parameters correspond to a bimodal Rb con- 
densate in \F — l,mp — 1) and \F = 2,mF = —1) with 
Na = ./Vfc = 5 X 10^ and where a Feshbach resonance is 
used to reduce Uab by about 10% with respect to its bare 
value [l4[|15j . The considered harmonic trap is very steep 
w = 27r X 2 kHz. In the figure we compare our modulus- 
phase approach (dashed line) with the full numerical so- 
lution (solid line) and with a stationary calculation using 
(|27| - (|28|) (dash-dotted line) which is equivalent to a two- 
mode model. The oscillation of the squeezing parameter 
in the two dynamical calculations (dashed line and solid 
line) are due to the fact that the sudden change in the 
mean-field causes oscillations in the wave functions whose 
amplitude and the frequency are different for each Fock 
state. From the figure, we find that our modulus-phase 
approach obtained integrating 5 Gross-Pitaevskii equa- 
tions (dashed line) reproduces the main characteristics of 
the full numerical simulation using 3000 Fock states (solid 
line). The stationary two mode model on the other hand is 
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- Two-mode stationary 

- Modulus-phase 5 GPE 

- Full numerical simulation 




2.5 3 
t[ms] 

Fig. 2. Spin squeezing as a function of time. Comparison of the 
modulus-phase model (red dashed line) with a full numerical 
calculation with 3000 Fock states (blue solid line) and with a 
stationary two-mode model (violet dash-dotted line). Spatial 
dynamics is strongly excited after the 7r/2-pulse populating a 
second internal state, lo = 2tv x 2 kHz, Na — Nb — 5 x 10*, 
m=87 a.m.u., a^a = 100.44rB, ati, = 95.47 rs, Uab ~ 88.28r_B. 
No particle losses, rs is the Bohr radius. 



not a good approximation in this case. Only for some par- 
ticular times the three curves almost touch. At these times 
the wave functions of all the Fock states almost overlap 
and, as we will show in our analytical treatment, spatial 
dynamics and spin dynamics disentangle. 

In FigOwe move to a shallow trap and less atoms. We 
note that in this case both the modulus-phase curve and 
the numerical simulation are very close to the stationary 
two-mode model which is then a good approximation at 
all times. 



Two-mode stationary 
Modulus-phase 5 GPE 
Full numerical simulation 




0.2 , , ,0.3 

t[s] 



X7l 



s- 




Fig. 3. (Top) spin squeezing as a function of time in a case in 
which the spatial dynamics is weakly excited. Blue solid line: 
full numerical calculation with 1000 Fock states. Red dashed 
line: modulus-phase model. Violet dash-dotted line: station- 
ary two-mode model. (Bottom) angle giving the direction of 
the mean spin projection in the equatorial plane of the Bloch 
sphere. Parameters: tj = 2-7r x 42.6 Hz, iVa = iV;, = 1 x lO", 
m=87 a.m.u., Uaa = 100.44r-B, att = 95.47 r_g, flab ~ 88.28rs. 
No particle losses. r_g is the Bohr radius. 



with 



2.4 Squeezing in the breathe-together solution 

In this section we restrict to a spherically symmetric har- 
monic potential JJ = vnuP'r^ (2 identical for the two 
internal sates. For values of the inter particle scattering 
lengths such that 

a-ab < aaa,abb (30) 

and for a particular choice of the mixing angle such that 
the mean field seen by the two condensates with Na and 
Nb particles is the same: 

Nagaa + NhQab = Nbgbb + NaQab = Ng , (31) 

the wave functions (pa and (pi, solve the same Gross-Pitaevskii 
equation. In the Thomas- Fermi limit, the wave functions 
(pa and (pb share the same scaling solution (p p^SllTO] and 
"breathe-together" [9]. 



(PC _ 



M 



c^h 



£4 

15 



-uj^L: 



1/2 r 



r2 1 ^/^ 

^0- 



(33) 

(34) 

(35) 



{r.t) 



iriit) 



C^/^{t) 



e™'''^(*)/2'^^(*)^o(r//:(t)) 
(32) 



fi is the chemical potential of the stationary condensate 
before the tt/2 pulse, when all the TV atoms are in state a, 
and Rq = y^2p,/muj'^ is the corresponding Thomas-Fermi 
radius. The initial conditions for (p4|) are C{0) = 1 and 
£(0) = 0. 

Note that the scaling solution identical for the two 
modes a and b is valid only for Na — Na, Nf, — Nb and 
does not apply to all the wave functions <paiNa, Nb) and 
(l)b{Na,Nb) in the expansion equation ^. Nevertheless, 
an advantage of choosing the mixing angle in order to sat- 
isfy the breathe-together condition equation (l3T|) . is that 
the mean spin has no drift velocity. In Fig|3]we calculate 
the spin squeezing (Top) and the angle if giving the direc- 
tion of the mean spin projection on the equatorial plane of 
the Bloch sphere (Bottom), for the same parameters as in 
FigEl except for the mixing angle that we now choose sat- 
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isfying equation ([3T|l while in Fig [3] we had Na — Nf,. Note 
that if practically does not evolve. The maximum amount 
of squeezing is lower in the breathe-together configuration 
than in the even-mixing case (see also [H]). However, as 
we will see in the next section, this conclusion does not 
hold when particle losses are taken into account. 
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Fig. 4. (Top) spin squeezing in breathe-together conditions as 
a function of time. Blue solid line: full numerical calculation. 
Red dashed line: modulus-phase model. Violet dash-dotted 
line: stationary two-mode model. (Bottom) angle giving the 
direction of the mean spin projection on the equatorial plane 
of the Bloch sphere. Parameters: uj ^ 27vx 42.6 Hz, Na = 7432, 
Nt = 12568, m=87 a.m.u., aaa = 100.44rs, abb = 95.47 rs, 
ciab = 88.28 rB. No particle losses, rs is the Bohr radius. 



By linearization of <j)a{Na, Nj,) and (pbiNa, Ni,) around 
the breathe-together solution (f> and using classical hydro- 
dynamics, it is even possible to calculate analytically the 
parameters Xd and Xs relevant for the squeezing dynamics 
[9]. One obtains: 



Xdir,t) 



2h\5N , 



9bb — '^9ab 



9a 




dt' 5 lmB{t) 

£3(?)^2 n^ 



Xs{r,t)^{\Cb\'-\Ca\')xd^Xo{r,t). (37) 



with 



^5- 



NaNb 9aa + 9bb - '^gab 



iV2 



9aa 



1/2 



5'/\ 



(38) 



and where ImS(t) is solution of the differential equations 

(39) 
(40) 






cHty 

to be solved together with equation (|34|) . with initial con- 
ditions ^(0) = B{0) = 1. In practice, when we expand the 
condensate wave functions around the breathe-together 
solution equation p2p as in [5], we encounter the hydro- 
dynamics operator S [20] 



SI 



_ Nga 



M 



'-div[(f>Q grad a] 



(41) 



The deviation of the relative phase and the relative den- 
sity from the breathe-together solution expand over two 
eigenmodes of S: A zero-energy mode which grows lin- 
early in time and gives the dominant features of phase 
dynamics and squeezing (integral term in the curly brack- 
ets in Eg. ([551) ). and a breathing mode of frequency H^ 
which is responsible for the oscillations of the squeezing 
parameter. The fact that in breathe-together conditions 
and within the modulus-phase approximation Xs = Xo is 
shown in Appendix [Cl 

We give an example corresponding to strongly oscil- 
lating wave functions in Fig[5] where we compare the spin 
squeezing from the analytical theory with a numerical 
simulation. In the analytical formula, the entanglement 
between spatial degrees of freedom and spin dynamics is 
apparent as Xd equation (j36p is position dependent. The 
points in which the dynamical curve (dotted line) touches 
the stationary two mode curve (dash-dotted line) corre- 
spond to ImB = (see the bottom curve) where space 
and spin dynamics are disentangled. We note however 
that the validity conditions of classical hydrodynamics are 
more stringent for a mixture of condensates with rather 
close scattering lengths than for a single condensate [S]. 
We checked numerically that in order for equation (j36p 
to correctly predict the frequency of the oscillations in 
the squeezing parameter, we have to enter deeply in the 
Thomas-Fermi regime. 



2.5 "Extracted" spin squeezing 

As we pointed out, the definitions equations pHll - p^ ex- 
plicitly include the spatial overlap between the two modes. 
Here we give an alternative definition that can be used al- 
ways, whether or not the modes overlap. To this aim, we 
introduce the time-dependent operators 



a= I d'rrair,t)Mr), 
b^ I d^r4>l{r,t)Mr), 



(42) 
(43) 



where 0e(r, t) is the solution of Gross-Pitaevskii equation 
(I7|) for mode e with Na, Nb particles. We then introduce 
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Fig. 5. (Top) test of the analytical formula equation (|36p in 
the deep Thomas-Fermi regime. Spin squeezing as a function 
of time. Blue solid line: full numerical calculation. Red dashed 
line: modulus-phase model. Black dotted line: analytical curve 
using equation (|36p . Violet dash-dotted line: stationary two- 
mode model using (I27|l - (|28p . (Bottom) function ImS(i). Spatial 
and spin dynamics disentangle when ImS(i) = 0. Parameters: 
Na = Nb = 5x 10^, uj = 2tt x2 kHz, m—87 a.m.u., Uaa = att = 
0.3 aho, lab = 0.24 a^o. a^o is the harmonic oscillator length: 
No particle losses. 
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the spin operators: 






(44) 
(45) 
(46) 



In the new definition of spin squeezing calculated by the 
spin operators defined in equations (|ii)) - ([i^ . which we 
call the "extracted" spin squeezing, we still take into ac- 
count entanglement between external motion and spin dy- 
namics, but we give up the information about the over- 
lap between the two modes. In Appendix [Dl we give the 
quantum averages useful to calculate the extracted spin 
squeezing within the modulus-phase approach described 
in Section 12.31 We will use this extracted spin squeezing 
in Section [5l 

Comparing the expressions given in Appendix |D] with 
those of Appendix iFl (in the absence of losses), one real- 
izes that in the stationary case, where Xd: Xs and xo £^re 
space independent, the extracted spin squeezing dynami- 
cal model reduces to a two-mode model that we study in 
detail in the next section. 



3 Two-mode model with Particle losses 

In this section we generalize our results of [12] to possibly 
overlapping and non-symmetric condensates. In subsec- 
tion EH] we address the general case, while in subsection 
13. 21 we restrict to symmetric condensates and perform an- 
alytically an optimization of the squeezing with respect 
to the trap frequency and number of atoms. In the whole 
section, as in [12j . we will limit to a two- mode stationary 
model and we do not address dynamical issues. 



3.1 Spin squeezing in presence of losses 

We consider a two-component Bose-Einstein condensate 
initially prepared in a phase state, that is with well defined 
relative phase between the two components, 

l^(0»Hri."^°l''"''°'Xi^'l'"''^')"|0)^ ,47, 



im 



When expanded over Fock states, the state pT]) shows bi- 
nomial coefficients which, for large A^, are peaked around 
the average number of particles in a and 6, Na and Nh- In 
the same spirit as the "modulus-phase" approximation of 
subsection 12. 31 we can use this fact to expand the Hamil- 
tonian of the system to the second order around Na and 
Afc 

i/o - S(Aa, Afc) + ^ y.,{N, - N,) + ^^N^^i,{N, - N,f 

+ I [dN.t^a + 9jV„/^b) {Na - Na){Nb - Nb) (48) 

where the chemical potentials /^^ and all the derivatives 
of ^e should be evaluated in Na and Nb- We can write 

H^ = ffj + hvf,{Na-Nb) + %Na-Nb? (49) 



with 



X 



1 

2h 
1 

2h 
1 

2h 



[tia - Mb) - fix{Na - Nb) + hx{N - iV)J(50) 

{dN^tJ-a + dNj-ib - dN^fla - ^Na^J■b) n^^n^ (51) 



(dNafJ-a - 9N^flb)N^,Nt 



(52) 



The function / of the total number of particles, N — 
Na + Nb, commutes with the density operator of the sys- 
tem and can be omitted. The second term in equation 
(j49)) proportional to Sz describes a rotation of the aver- 
age spin vector around the z axis with velocity Vj^. The 
third term proportional to S"^ provides the nonlinearity 
responsible for spin squeezing. It also provides a second 
contribution to the drift of the relative phase between the 
two condensates in the case Na ^ Nb- 

In presence of losses, the evolution is ruled by a master 
equation for the density operator p of the system. In the 
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interaction picture with respect to Hq, with one, two, and 
three-body losses, we have: 



3 



l = EE7i™^ 



dt 



m=l e=a,b 
lab 



cT~pc\-^-\{c^rcT,~p} 



CaCbPcicl 



TilciclcaCb, p} 



(53) 



where p = e^"°*/^pe-'"°*^'', Ca = e'"°^/''ae-'"'>*^''\ and 
similarly for b, 



v(™) 



K, 



(e) 



m 

Kab 



d'r\Mr)r, 



lab^^ I d'r\Urr\Mr)\' 



(54) 
(55) 



Km is the TO-body rate constant {m = 1,2,3) and 0e(r) 
is the condensate wave function for the e component with 
Na = Na and Nb = Nb particles. Kab is the rate constant 
for a two-body loss event in which two particles coming 
from different components are lost at once. 

In the Monte Carlo wave function approach [21] we de- 
fine an effective Hamiltonian iJcfr and the jump operators 



state and we can calculate exactly the evolution of the 
state vector analytically, as we did in J12j for symmetrical 
condensates. When two and three-body losses enter into 
play, we introduce a constant loss rate approximation |22j 



m=l e=a,b 

valid when a small fraction of particles is lost at the time 
at which the best squeezing is achieved. In this approxi- 
mation, the mean number of particles at time t is 



{N)=N 1- EE^i'^^+^^'^h (62) 

\s:=a,b m / 

ri") = 7V™-im7i") ; Fab = labV^Nb (63) 

where for example F^ t is the fraction of lost particles 
due to TO-body losses in the e condensate. Let us present 
the evolution of a single quantum trajectory: Within the 
constant loss rate approximation, we can move all the 
jump operators in ([58]) to the right. We obtain: 



\nt)) 



-At/2 



nJirHi.)i'^(o)) 



H, 



ih 



off - 



E E 7i"'ct^ 



ih 



J=l 



ni—l e—a.b 



Ji™^ = V7- 



,("0, 



7(2) , — 

Jab = Vl^CaCb 



labc\clcaCb (56) 



(57) 



->^''\-'^>'\ak\W + (3k)N-N[k) (64) 



where 



We assume that a small fraction of particles will be lost 
during the evolution so that we can consider x, le and 
lab as constant parameters of the model. The state evolu- 
tion in a single quantum trajectory is a sequence of ran- 
dom quantum jumps at times tj and non-unitary Hamil- 
tonian evolutions of duration tj : 

. . . 4™i)(ii)e-'-f^'=«^i/'"'|if (0)) , (58) 

where now Sj — a,b oi ah. Application of a jump J^j ' (tj ) 
to the A^-particlc phase state at tj yields 

j'i^'\h)\^)N ex 10 + A,t,)t,^ra, , (59) 

A- = '^X^ej,ab + (X + X)fnjS^.^a + {x ~ X)"lj ^e, ,6(60) 

After a quantum jump, the phase state is changed into a 
new phase state, with m particle less and with the relative 
phase between the two modes showing a random shift Ajtj 
with respect to the phase before the jump. Note that in 
the symmetrical case x = and no random phase shift 
occurs in the case of a jump of ah. Indeed we will find 
that at short times in the symmetrical case theses kind of 
crossed ab losses are harmless to the the squeezing. 

In presence of one-body losses only, also the effective 
Hamiltonian changes a phase state into another phase 



i"'=i' = n E E '-. 

j=l \^m' = 1,2,3 e'=a,b 

+ N^\Ca\^\CbflabSe,,ab} 
k 



,m'5e,,e'N"^'\C,.\^"''l'^f^ 



(65) 



i'5„ 



(3k^Y.^A2x5e,,ab+ Y. 

j = l [ m' = l,2,3 

X [(X + X)Se, ,a - (X - X)Se, ,b] } 
k 

^(^) ^Y. Y. ^'Srn„m'{Se,,a + ^e,,b) + 26,^, ab , (67) 
j = l m' = l,2,3 



(66) 



and Tk is a phase which cancels out when taking the av- 
erages of the observables [33] . 

The expectation value of any observable O is obtained 
by averaging over all possible stochastic realizations, that 
is all kinds, times and number of quantum jumps, each 
trajectory being weighted by its probability 

{0)^Y.I dtidt2---dtk J2 (•^WI'^I'^W)- (68) 

Note that the single trajectory ([M)) is not normalized. The 

prefactor will provide its correct "weight" in the average. 

We report in Appendix |E] and [F] the averages needed 

to calculate the spin squeezing for one-body losses only 
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(exact solution) and for one, two and three-body losses 
(constant loss rate approximation) respectively. The ana- 
lytical results are expressed in terms of the parameters x 
and X defined in equations (|5T|) and ([52|l respectively and 
of the drift velocity 

V - ^ [(Ma - Mb) - MNa - m) + hx{N - N)] , (69) 
where N is the total initial number of atoms. 



3.2 Symmetrical case: optimization of spin squeezing 

If we restrict to symmetrical condensates which may or 
may not overlap, we can carry out analytically the opti- 
mization of squeezing in presence of losses. In the symmet- 
ric case and constant loss rate approximation it turns out 
that ASl = (iV)/4. This allows to express f^ in a simple 
way: 



e 



(g^g) 
(6ta)2 



({a)a)+A-^/WT&^) , (70) 



with 



i= -Re{{h'^a^ah^b'^b^aa)) 
B = 21m{(b^b''ba)) . 



(71) 
(72) 



An analytical expression for spin squeezing is calculated 
from ([701) with 



(6ta) 



-At 



cos^^-^(xt)A^^i 



(73) 



-At 



A 



B = 



-N{N - 1) [Fa ~ F2 cos^-2(2^<)] (74) 



,-At 



■ COS 



W-2 



{Xt) sinixt)NiN - l)Fi (75) 



where the operator N — {N — da) acts on the functions 



Ff3{a) = exp 



3 



2j^"'he' 



, sin(m/3xO 
m(3xi cos™ {Pxt) 



+ 



labte 



la 



cos'^{Pxt) 



(76) 



with /3 = 0, 1, 2 and all expressions should be evaluated in 
(7 = In TV. 

We want now to find simple results for the best squeez- 
ing and the best squeezing time in the large N limit. In 
the absence of losses [1] the best squeezing and the best 
squeezing time in units of 1/x scale as N~^^^. We then 
set N ~ e^^ and rescale the time as xt = ''"S^- We expand 
(|70l) for e <C 1 up to order 2 included, keeping F^^^^x con- 
stant. The key point is that in this expansion, for large N 
and short times, the crossed losses ab do not contribute. 
As in [12], introducing the squeezing ^o(t) in the no-loss 
case, we obtain: 



eit)=eoit) 



1 r,„t 

1 + -- ^ 



Hiit) 



(77) 



ith: 



^sq = E^q"^^^dr4" 



mT^"''> 



(78) 



The result (|77|) very simply accesses the impact of losses 
on spin squeezing. First it shows that losses cannot be ne- 
glected as soon as the lost fraction of particles is of the 
order of ^q. Second it shows that in the limit N -^ 00 
and ^o(^bcst) — ^ 0, the squeezing in presence of losses is of 
the order of the lost fraction of particles at the best time: 
C^(^bcst) ~ ^sqibcst/3. This also sets the limits of validity 
of our constant loss rate approximation. For our approx- 
imation to be valid, the lost fraction of particle, hence 
squeezing parameter at the best squeezing time, should 
be small. 

From now on, the optimization of the squeezing in the 
large N limit proceeds much as in the case of spatially 
separated condensates [12] • The only difference is in the 
stationary wave functions in the Thomas- Fermi limit. For 
overlapping condensates we consider a stable mixture with 



aab < a-aa — abb : 



(79) 



and we introduce the sum and difference of the intra and 
inter-species s-wave scattering lengths: 



ad — aaa ~ aab ■ 



In the symmetric case considered here we have 



Ma 



X 



fib = -huj 



15Nas 



33/532/5 ^ f^ 
M 



53/5 



2 a/io 

1/5 



2/5 



-6/5^-3/5^ 
3/5 

as 



(80) 
(81) 



(82) 

(83) 
(84) 



J-(2) ^ 
7-(3) = 



-6/5 



2^(^] Cu'/^N^/^a;'/'K2 (85) 



54/5 



219/531/57^2 ypj 



-12/5 

— ] ^^12/5^4/5 



(86) 



where aho = y^h/AIuj is the harmonic oscillator length, u) 
is the geometric mean of the trap frequencies. We recover 
the case of spatially separated condensates [Hj setting 
Oab^Oin dHOl-dll]). 

The squeezing parameter for the best squeezing time 
^'^{thcst,^) is minimized for an optimized trap frequency 



ui 



,opt 



219/1275/12^5/6 ^ ^1/2 /^ X5/12 



151/3 M7V1/3 \K3 



(87) 



Note however that this optimization concerns one- and 
three-body losses only. The effect of decoherence due two 
two-body losses quantified by the ratio /"^^Vx i^ indepen- 
dent of the trap frequency. 
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Once the trap frequency is optimized, C^(ibcst,'^°^*) 
is a decreasing function of N. The lower bound for ^^, 
reached for iV ^ oo is then 



inf r 

t,u),N 



5V3A/ 
287rfi 



2/3 



2 aH Qd 



2/3 



A simple outcome of this analytic study is that, for pos- 
itive scattering lengths Uaa, a-ab, the maximum squeezing 
is obtained when aab = that is for example for spatially 
separated condensates. Another possibility is to use a Fes- 
chbach resonance to decrease the inter-species scattering 
length aab [H1IT5] , knowing that the crossed a — b losses 
do not harm the squeezing at short times. 



4 Results for overlapping condensates 

In this and the next section wc give practical examples of 
application of the analysis led in the two previous sections. 



4.1 Feshbach resonance-tuned bimodal Rb BEC 

We consider a bimodal Rb condensate in \F = 1, mp — 1) 
and \F = 2, nip = —1) states where the scattering length 
Qab is lowered by about 10% with respect to its bare value 
using a Feshbach resonance p!4l[T5] . 

In Figl6] (Top) and (Bottom) we compare a situation 
in which the initial condensate is split evenly in the a and 
b components to a situation in which the mixing is cho- 
sen in order to satisfy the "breathe-together" conditions 
(PT|) . For the considered parameters, which are the same 
as Figl3] and FiglH the spatial dynamics is not important 
and the two-mode model is a good approximation at all 
times. 

The squeezing in presence of losses is calculated using 
our general results of Section 13.11 for asymmetric conden- 
sates. Although without losses the even splitting is more 
favorable, with one, two, three-body losses results are com- 
parable ^^ ~ 6 X 10~^. We also show a curve obtained for 
one and three-body losses only (dashed-line). It is clear 
that for the considered Rb states the dominant contribu- 
tion for decoherence comes from the two-body losses in 
the F ~ 2 state severely limiting the maximum amount 
of obtainable squeezing. 

In the cases considered in Fig|S] asymmetric two-body 
losses are very high, we therefore check the validity of 
the constant loss rate approximation with an exact Monte 
Carlo wave function simulation. The main result is that 
the constant loss rate approximation is accurate up to 
the best squeezing time. A more complete discussion is 
presented in Appendix [Gl 
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Fig. 6. Spin squeezing with and without losses in a bimodal 
Rb condensate from the stationary two-mode model. A Fesh- 
bach resonance is used to reduce the inter-species scattering 
length. Violet dash-dotted line: without losses. Blue dashed 
line: with one and three-body losses. Red solid line: with 
one, two and three-body losses. (Top) with a 50 — 50 mix- 
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of using Na atoms in the \F — l,mp — ±1) states [5]. 
In Figlllwe calculate the best obtainable squeezing with 
these two states. Parameters are chosen according to our 
optimization procedure of Section 13.21 A large amount of 
squeezing ^^ — 1.9 x 10^"^ can be reached at the best 
squeezing time. 



4.2 Bimodal BEC of Na atoms 

By using two states in the lower hyperfine manifold, one 
can greatly reduce two body losses. A possible example is 



Using our full numerical and our approximated dynam- 
ical approaches, (not shown) we checked that the two- 
mode model is an excellent approximation for these pa- 
rameters. 
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Fig. 7. Spin squeezing with and without losses in a bimodal Na 
condensate from the stationary two-mode model. Violet dash- 
dotted line: without losses. Blue dashed line: with one and 
three-body losses. Optimized parameters: Na = iV;, = 4 x 10* 
(jj — 2tt X 183 Hz, m—23 a.m.u., aaa = att = 51.89 rs, Uat ~ 
48.25 tb, rs is the Bohr radius, x = 5.517 x 10"^s~\ X = 0. 
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5 Dynamically separated Rb BEC 

In this subsection we consider a bimodal Rb condensate 
in \F = l,mp = —1) and \F = 2, nip = 1) states. Rather 
than using a Feshbach resonance to change gab, we con- 
sider the possibility of suddenly separating the two clouds 
right after the mixing 7r/2 pulse using state-selective po- 
tentials [28], and recombining them after a well chosen 
interaction time. A related scheme using Bragg pulses in 
the frame of atom interferometry was proposed in [29] . We 
consider disc shaped identical traps for the two states a 
and b with ujz > ujx,y = ^±, that can be displaced inde- 
pendently along the z axes. In order to minimize center- 
of-mass excitation of the cloud, we use a triangular ramp 
for the displacement velocity, as shown in Fig[8](Bottom), 
with total move-out time 2r = 4:n/iVz |30]. In Fig|8](Top) 
we show the z-dependence of densities of the clouds, in- 
tegrated in the perpendicular xy plane, as the clouds are 
separated and put back together after a given interaction 
time. 

We use our dynamical modulus-phase model in 3 di- 
mensions to calculate the spin squeezing in this scheme. 
As the spatial overlap between the two clouds reduces 
a lot as they are separated, in Fig|9] we calculate both 
the spin squeezing obtained from the definitions ([T0| - (fT2|l 
of spin operators (dashed line), and the "extracted spin 
squeezing" introduced in Section[23]based on the "instan- 
taneous modes" (|^^ - P5)) (solid line). The oscillations in 
the dashed line are due to tiny residual center of mass os- 
cillations of the clouds that change periodically the small 
overlap between the two modes. They are absent in the 
extracted spin squeezing curve (solid line) as they do not 
affect the spin dynamics. When the clouds are put back 
together and the overlap between the modes is large again, 
the spin squeezing and the extracted spin squeezing curves 
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Fig. 8. (Top) \cj>a{z,t)f and [0f,(2;,i)|^ in arbitrary units as 
the clouds are separated and put back together after an inter- 
action time of about 15 ms. The harmonic potential for the 
a-component does not move, while that for the 6-component 
is shifted vertically with a speed 5z. The distance between the 
two trap centers when they are separated is 5z = Ay/h/Mui^. 
(Bottom) variation in time of &z. Parameters: Na = iV(, = 
5 X IC*, Ld^^y = 27r X 2.31 Hz, ci;^ = 27r x 1 kHz, m=87 a.m.u., 
daa ~ 100.44 rs, abb ~ 95.47 rs, aab = 98.09 rs, rs is the Bohr 
radius. No particle losses. 



give close results (not identical as the overlap of the two 
clouds is not precisely one). 

In Fig llOl (Top) we compare the extracted spin squeez- 
ing curve of Fig[9] (solid line) with a two-mode stationary 
calculation (dash-dotted line) assuming stationary con- 
densates in separated wells. We notice that the squeezing 
progresses much more slowly in the dynamical case. In- 
deed when we separate the clouds, the mean field changes 
suddenly for each component exciting a breathing mode 
whose amplitude and frequency is different for each of the 
Fock states involved. In the quasi 2D configuration con- 
sidered here, the breathing of the wave functions is well 
described by a scaling solution in 2D for each condensate 
separately [TSlfTO] adapted to the case in which the trap 
frequency is not changed, but the mean-field is changed 
suddenly after separating the two internal states: 



'Xr^,t) = 






^imriC,(,t)/2hC,{t)^^ 



Ce{t) 



(89) 
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Fig. 9. Spin squeezing as the two Rb condensates are separated 
and put back together after an interaction time of about 15 ms. 
Red dashed line: Spin squeezing obtained from the definitions 
(|10p - (|12p of the spin operators explicitly including the overlap 
between the clouds. Blue solid line: extracted spin squeezing 
based of the "instantaneous modes" (|42|l - (|43[l . Parameters as 
in FigEl 
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(91) 

(92) 



H is the chemical potential of the stationary condensate 
before the tt/2 pulse, when all the N atoms are in state 
a, i?o — \/2pJrnuj^ is the corresponding Thomas-Fermi 
radius, and g^^ is a reduced coupling constant to describe 
the interaction between two atoms in the e condensate in 
quasi 2D system, where we assume that the condensate 
wave functions in the confined direction are Gaussians: 



ge 






27rS 



(93) 



with flej the 3D scattering length. The initial conditions 
for dSI]) are £^(0) = 1 and £^(0) = 0. 

We can use ([M|) to calculate the squeezing (dotted 
curve) and we note that it reproduces well the spin squeez- 
ing curve obtained integrating 5 Gross-Pitaevskii equa- 
tions in 3D (full line). As we studied in detail in Section 
12.41 oscillations of the wave functions cause oscillations 
of the squeezing parameter due to entanglement between 
spatial and spin dynamics. Indeed what we see in the ex- 
tracted spin squeezing curve of Fig[TO] (Top) is the begin- 
ning of a slow oscillation for the squeezing parameter. In 
FigdO] (Bottom) we show the long time behavior. There 
are indeed times at which the spatial and spin dynamics 
disentangle, and the dynamical curve and the steady state 
curve touch (see Sect. 12.4) ). Unfortunately these times are 



plO" 



10 



Extracted spin squeezing 
Scaling solution 
Stationary solution 
Stationary solution with 
1 ,2,3-body losses 




2 4 6 8, JO, 12 14 16 1{ 
t [ms] 



10" 



10" 



nIO" 



10" 



10" 





1 1 1 


' 






— Scaling solution 






Stationary solution 








Stationary solution witin 








1,2,3-body losses 




_________ — ■- 


\ 
\ 


^^ 


\ '. .-■■ ■--. 


.-.'--' 


^ ■, .' ■■--. 


---" 


\ \ .' ."' 




"^ ■ ' ,' _, . -^ 




X \ • ^ . -' 


■ 















0.05 



0.1 



0.15 
t[s] 



0.2 0.25 



Fig. 10. Spin squeezing as a function of time. (Top) compari- 
son between a dynamical calculation and a stationary calcula- 
tion. Blue solid line: extracted spin squeezing in 3D. Black 
doted line: 2D scaling solution based on (|9ip . Violet dash- 
dotted line: stationary calculation in 3D without losses. Red 
solid line: stationary calculation in 3D with losses. Spin squeez- 
ing progresses more slowly in the dynamical calculation than 
in the stationary calculation. (Bottom) long time behavior. 
Black doted line: scaling solution. Violet dash-dotted line: sta- 
tionary calculation without losses. Red solid line: stationary 
calculation with losses. Parameters: x ~ 5.003 x 10~^s~^, 
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The other parameters are as in Fig|8] 



not accessible here in presence of losses (in particular the 
high two-body losses in the higher hyperfine state) . Notice 
that in the first 15 ms of evolution considered in Fig|H]and 
FigdO] (Top) the effect of losses is small and the main lim- 
itation at short times is provided by the spatial dynamics. 

For a lower number of atoms, the sudden change in the 
mean field and the consequent oscillations of the squeez- 
ing parameter are reduced. In Fig llll we show the spin 
squeezing obtained by suddenly separating two BEG of Rb 
atoms in \F = l,mF = —1) and \F = 2,mF = 1) states 
with 1000 atoms in each component. The dotted line is 
a dynamical calculation using the quasi 2D scaling solu- 
tion (pij) (and no losses), while the dash-dotted line and 
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Fig. 11. Spin squeezing as a function of time in two smaU 
Rb condensates. Black doted line: scaling solution based on 
l|9ip . Violet dash-dotted line: stationary calculation without 
losses. Red solid line: stationary calculation with losses. Pa- 



(a) 



K\ 



rameters: K\ 
10-2im3s-i[24],A-f' 



(6) 



,(a) 
-1 



Ab) 



0.01s"\ K^2 

The other param- 
eters: Na = Ni = 10^, uj^^y = 27r X 11.82 Hz, oj^ = 27r x 2 kHz, 
m=87 a.m.u., aaa — 100.44rs, abb = 95.47 rs, aab = 98.09r-s, 
rs is the Bohr radius, x = 0.213s"\ x = 2.763 x 10"^s"\ 



with M. Oberthaler, J. Esteve and K. M0lmer. Our group 
is a member of IFRAF. 



A Quantum averages of the field operators 



Using equations (l20|)-([2T]l. the averages needed to calcu- 
late squeezing parameter can be written in terms of the 
wave functions (/)a, 06 and the phase factor A solution of 
equation Q: 



-ACal'^^^-'^lar^C^Ca 



xc^l{Na-l,Nb+l, r)^a {Na,Nb,r) 

X exp{i[AiNa - 1, iV5 + 1) - A{Na,Nt)]/h} 

X [{MNa -l,Nb + l)\MNa, iVb))]^»-l 

X [{MNa - 1, iVb + l)\MNa, Nb))]''' . (94) 



the solid line are stationary calculations without and with 
losses respectively. Note that around t — 0.02s, where the 
dynamical curve and the stationary curve touch, a squeez- 
ing of about ^^ ~ 2 X 10^^ could be reached despite the 
high losses in the F = 2 state |31] . 



6 Conclusions 

In conclusion we developed a method to study the en- 
tangled spatial and spin dynamics in binary mixtures of 
Bose-Einstein condensates. The method, which is the nat- 
ural extension of our work [9] to the case of spin squeezing, 
allows a full analytical treatment in some cases and can 
be used in the general case to study a priori complicated 
situations in 3D without the need of heavy numerics. In- 
cluding the effect of particle losses and spatial dynamics, 
we have calculated the maximum squeezing obtainable in 
a bimodal condensate of Na atoms in |_F = 1,tof = ±1) 
states when the two condensates overlap in space, and we 
have calculated the squeezing in a bimodal Rb conden- 
sate in which a Feshbach resonance is used to reduce the 
inter-species scattering length as recently realized experi- 
mentally [15]. For Rb we also propose an original scheme 
in which the two components are spatially separated us- 
ing state-dependent potentials, recently realized for the 
\F = l,mF = —1) and \F = 2,tof = 1) states, and then 
recombined after a well chosen squeezing time. With this 
method we show that ^^ ~ 2 x 10~^ could be reached 
in condensates of 1000 atoms, despite the high two-body 
losses in the higher hyperfine state. 

Yun Li acknowledges support from the ENS-ECNU 
program, and A.S. acknowledges stimulating discussions 



E 



iV! 



nr^iCbn 



X(j)l{Na. Nb, r)^l{Na, Nb, r>a(iV„ Nb, r) 
xMNa,Nb,r'). (95) 



{^lir)i'lir')Mr)Mr')) 



Nl 



iNa-2y.Nb\ 



|^^|2(Ar.-2)|^^|2Ar,^*2^2 



N 

E 

Na=2 

xcj)l (Na -2,Nb + 2, r)rb{Na -2,Nb + 2, r') 

X(j)a{Na,Nb,r)(t)aiNa,Nb,r') 

X exp{z[A(iV, - 2, iVb + 2) - A{Na,Nb)]/h} 

X [{MNa -2,Nb + 2)\(baiNa, Nb))f^'^ 

X [{MNa -2,Nb + 2)[MNa,Nb))f' . (96) 



{i^l{r)^jl{r')Mr)Mr')) 

Nl 



"bV IVb 
N-l 



{Na - mm - ly. 



iap(^"-i)|C6p^^c:c, 



E 

xMNa -l,Nb + 1, r)cl)l{Na -l,Nb + 1, r') 

X(t)b{Na,Nb,r)MNa,Nb,r') 

X exp{i[A{Na ~l,Nb + l)- A{Na, Nb)]/h} 

x[{MNa-l,Nb + lMa{Na,Nb))f--^ 

X [{MNa -l,Nb + l)[MNa, Nb))f'-' . (97) 
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xMNa,Ni„r)cf>t{Na,m,r') 

X exp{i[^(Ar, + 1, TVb - 1) - A(iV„ iVb)]/n} 

X HMNa + l,m- l)\MNa, Nb))f'-^ . 



(98) 



We use these averages to calculate the squeezing in our full 
dynamical model. In practice we do not sum over all the 
Fock states but over a "large enough width" (typically 
> 6viV) around the average number of atoms Na, Nb. 
The spin squeezing is obtained by equation (fT5|) using the 
definitions p0l) - p2)) for the spin operators. 



B Quantum averages in the modulus-phase 
approach 



Within the modulus-phase approximation, the scalar prod- 
uct of the wave vectors can be written as 



{cl>a{Na -f3,Nb + (3)\MNa,Nb)) 

^ exp{if3 J d^r\4>a{r)\^[XQ{r) + Xd{r)]} 

{MNa-P,Nb+P)\MNa,Nb)) 

= exp{z/3/ dV|^,(r)nxo(r) - xd(r)]} 



(99) 
(100) 



{cl>b{Na-f3,Nb + /3)\MNa,Nb)) 
= J (fr4)l{r)^a{r) exp[i(iVa - I3)xd{r) - iNbXd{r)] 
X eMi{N - N)xs{r) - iNi\Ca\^ - \Cb\')xdir)] 
X exp[iPxo{r)] (101) 

{MNa+P,Nb~(3)\MNa,Nb)) 

= J d^rra{r)Mr) exp[-iNaXd{r) + i{Nb - P)xd{r)] 
X exp[-i(iV - N)xs{r) - iN{\Ca\^ - \Cb\^)xd{r)] 
xexp[-i/?xo(r)] (102) 

where /3 G Z, and we have used the relation 

J d^r\4>,\'' exp[i{dN^ - dN,)eeiNa,Nb)] 

c^exp[iJd^r\^,\^dN^-dN,)9,iNa,Nb)].il03) 

By using the Gross-Pitaevskii equations ([7]) for (pd^a, ^b) 
and for (j)^{Na, Nb), one obtains 



ihdt 



(Na - Na) 



86, 



+ [Nb ~ Nb) 



dNb 



Na,Nb 



dNa 

N,)g,em^ + [Ne' - N,.)gee'\<f>e'f , (104) 



where e y^ e' = a,b. Using (|104p together with the initial 
condition ([H]), we obtain for the phase factor A in Eq. ^ 



[A{Na -l,Nb + l)- A{Na, Nb)]/n 

HNa-l)Jd^r\Mr)\'[xo{r)+Xd{r)] 

-NbJd^r\Mr)\^[Xo{r)-Xd{r)] 



(105) 



[A{Na -2,Nb + 2)- A{Na, Nb)]/h 

= -2{Na - 2) / dMMr)nXo{r) + Xd{r)] 

-2Nbjd^r\4>b(r)\'[xdr)-Xd{r)] 

-jdM\Ur)nx^{r) + Xd{r)] 

+ \4>b{r)\''[xdr) - xdir)]} 

The averages and variances of the spin operators equa- 
tions (Unil-dlll) are obtained by equations (HHl-dMl) after 
spatial integration. We get: 

j d^r{i:l{r)Mr)) 
= NC*bCaSd^r4>l{r)4>a{r)[\Ca\'e'^-^-^ + |Cb|2e-*x.W]JV-i 



(106) 



X exp[i(7V - N)xs{r)] exp[-iiV(|CaP ■ 
xexp[Jxo(r-)]. 



n?)Xd{r)] 



(107) 



i d'^rd^r' {i^l{r)i,l{r')i,a{r)Mr')) 

^N{N-l)\Ca\'\Cb\'jd\dy4'l{r)Mr)ra{r')Ur'). 

(108) 

jd^rd^r'{^l{r)^l{r')^a{r)Mr')) 

= N{N~l)C*b^Clid^rd^r'4>l{r)4>a{r)rb{r')Ur') 

X\\C |2g«Xd('')+«Xd(i-') ^ \(j ]2^-iXd(r)-iXd(r'hN-2 

xeM-^N{\Ca\'~\Cb\')[xd{r)+Xd{r')]}eMnXo{r) 

+Xo{r')]}eMi{N~N)[xs{r)+Xs{r')]} 

X exp{-i / dV'(|0,|2[xo(r") + X<;(r")] 

+ m''[xo{r") - Xd{r")])) . (109) 

jd^rd^r'{^l{r)^l{r')Mr)Mr')) 
= N{N -l)C*bCa\Cb\^ j d^r'4>*b{r')Ur')[\Ca?e^^-^-') 
+ \Cb\^e-^^-'^^">f-^eM-^N{\Ca\^ - \Cb\^)Xd{r')] 
X exp[zxo(r-') - iXd{r')] exp[z(A^ - N)xs{r')] ■ (110) 

jd^rd^r'{'^l{r)i,l{r')^a{r)Mr')) 
= N{N -l)C:Cb\Ca? j d^r'4>l{r')Mr')[\Ca?e'^'''^-') 

MCb^'e^^-^^'y-^ eM^N{\Ca\' ~ \Cb\')Xd{r')] 

X exp[-ixo(r') - iXd{r')] exph^(Ar - N)xs{.r')] . (Ill) 

In the above expressions Xd, Xs and xo are the space and 
time dependent functions defined in equations p4| . (p5)) 
and (PS)1 . In practice it is sufficient to evolve five wave 
functions <j)a{r,t), (j)b{r,t) for {Na,Nb ± 6Nb) and {Na ± 
SNa, Nb) with 6Na,b 7^ (to calculate numerically Xd Xs 
Xo), and with SNa^b = (to calculate the central wave 
functions 4>a.b)- The spin squeezing is obtained by equation 
([T5|) using the definitions ([T0| - (fT2|l for the spin operators. 
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C Equality of Xs aid xo '" ^^^ 
breathe-together configuration 



Evaluating ((T04| ior e ^ a, Na ^ Na, e ^ b, Nb 
subtracting the two relations, on obtains 



d~9a d\ 



dN. dN„ 



= 



Nh and 



(112) 



where we used the fact that in breathe-together conditions 
1 0a I — 1 06 1- Equation (|112p implies that the time deriva- 
tive of Xs — Xo is zero. As for i = Xs = Xo = 0, we 
conclude that Xs = Xo at all times. 



D Extracted spin squeezing quantum 
averages 



{Wb-a) = N{N -l)C;Ca\Ctf J d'nd'r2\Mri)\^\Mr2)\^ 

X exp[ixS'^(r2,ri) - ixTir2,ri)] 
xexp[-tN{\Ca\'-\Cb\^)xTir2,r,)] 
xexp[t{N-N)xT{r2,ri)] (119) 

{a^d^~ab) ^ N{N -l)C:Cb\Ca\' J d^nd^r2\Mri)\^\Mr2)\' 

X exp[-ixS''(r-i, ra) - Jx^l^i- ^2)] 

xexp[iN{\Ca\'-\af)xT{ri,r2)] 

xexp[-i{N-N)xTiri,r2)] (120) 

In case the wave functions (pa-, </>& are stationary we re- 
cover the stationary two-mode model averages given in 
the next appendix in the particular case of no losses. The 
spin squeezing is obtained by equation (J15p using the def- 
initions (|44p - (P5|) for the spin operators. 



E Quantum averages with one-body losses: 

By using the instantaneous modes (I42l)-(l43]) and within Exact solution in the non symmetric Case 

the modulus-phase approach, the quantum averages useful ^ i- t • 1 

to calculate spin squeezing are expressed in terms of the ^^ ^^^^ appendix we give the exact result for quantum 
functions- averages needed to calculate spin squeezing in the case 

of a two-mode model with one-body losses only, in the 
general non-symmetric case. 



xT{r,r') = lid^^ - dNjMr) - 9b{r')]{N,,m) (113) 
xTir,r') = ^{dN^+dNj[ea{r)-eb{r')]{Na,m) (114) 



We obtain: 

{b^~a) = NC;Ca J d^nd^r2\Mri)\'\Mr2)\^ 

X exp[i(7V - iV)xr(r-2, ri) -I- ixS"(r-2, ri)] 
xexp[-iN{\Ca\'~\a\^)xT{r2,ri)] 



(atfltaa) = \Ca\''N{N - I) eM~^lat) 
(b^b^bb) = \Cb\^N{N-l)exp{~2-fbt) 



(121) 
(122) 
(123) 



{b^a^ab) = \Ct,\^\Ca\^N{N-l)exp[-{ja+Jb)t] (124) 



(b^d^ab)^N{N-l)\Caf\Cb\^ 

(b^'b^aa) = N{N - 1)C;^C^ d^rid^r2d^r3d^r4 
x\Mri)\'\Mr2r\Ur3r\Mr4)\' 

X exp{2i[xS''(r4,r2) +xS''(''3,r-i)]} 



(116) 
(117) 



(b^a) =C;Cbe-^'''*Nexp 

-2ivt 
1 



-l^{.la+lb)t 

(b^b^ba) = \Cbfc;Cae-^'^'N{N - l)e'''' 



-N-l, 



X exp 



-{"fa +3jb)t 



-N~2 



{o)a^ab) = |Cj2QCbe2™*iV(iV - l)e'x* 



X exp 



-l^{'ila+lb)t 



-N-2 



{b'^bh 



/~<*2/~<2 ^~4ivt 
^b ^a^ 



N{N ~1) 



X exp [-(7a -f 7fc)t]L^ ^ 



where we introduced the function Lp with /3 

\Ca? 



(125) 

(126) 

(127) 

(128) 
4,1,2 



L, 



X exp{-ziV(|C„|2 - \Cb\'')[xT{rA.r2)+xT{r,,r,)]} "^P " 7, + z/3(x + x) 

xexp{*(7V-iV)[xr(r-4,r2)+xr(r3,r-i)]} [gp 

xexp{-i/d3r5(|</.a(r5)|2[xg'^(r5,r5)-|-xS''('^5,r5)] ^ ^^ _ ^^(^ _ ^) 
+ \4>b{r,mxT{r,,r,)~x7{r,,r,)])} (118) 



7a e 



*''**-fi/3(x + x)e"^''"+'^''^* 



76e*''**-i^(x-x)e"''"'"''^^^* 



(129) 
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and V given by ([69 



F Quantum averages with one, two, 
three-body losses in the non-symmetric case 



( 3 



exp ■ 



E 



7^(1 



-im/3(x-x)tl 






imPix + x)[|CfcPe~''3(x-x)t + \Ca\^e'f>(x+x)t] 



oO-i+cr2 



lab[l 



^i2Pxt] 



z2/?x[|C6|2e-*'3(x-x)t + |C„|2gZ/3(x+x)t 



(139) 



In this appendix we give the quantum averages useful to 

calculate spin squeezing for the two-mode model in the 

general non-symmetric case, in presence of one, two and with /3 = 0, 1,2, and all the expressions should be eval- 
uated in a I — IniVfj, (72 = InA^f,. The expression of v is 
given in ()69|) . The spin squeezing is obtained by equation 
QoQ-i p^ using the definitions (Hil) - (|^ for the spin operators 
(with a = a and b = b). 



three-body losses. 



(atfltaa) = \Ca\^e-^' [N - {d^, + d„, - 1)] 
{b^bhb) = ICfcl^e-^* [N - (9,, + d„, - 1)] 



(131) 



G Test of the constant loss rate 
approximation for high asymmetric losses 



(132) 



{b'^a)ab) = ICfcHCape-^* [N - {d^, + d,, - 1)] 
X [N-{d„,+d„,)]Fa{ai,a2) 



X [N~{d,,+d,.,)]Fi{<Ji,a2) 



The constant loss rate approximation (I6ip is in general 
valid when a small fraction of particles is lost. In the case 
(133)of symmetric condensates, from equation (|77p one sees 
that the best squeezing in presence of losses is of the order 

(6^ a) = C*Cf,e^^2™+'*'^* i\C Pe"*'^* -f \Ci Pe*^*)^~^ *^^ ^^^ ^°^* fraction. So that ^(tbcst) ^ 1 guarantees that 

, „ . sthe lost fraction is small and the constant loss rate approx- 
^ imation is accurate. In the case of asymmetric condensates 
and asymmetric losses there might be other effects to con- 
sider as the population ratio between the two spin compo- 
nents might change in reality while it remains constant in 
the constant loss rate approximation. Indeed with the ap- 
X Fi(cri,o'2) (135)pi-oximation (|6ip. the initial phase state remains a phase 

state through out the whole evolution. As a consequence, 
_ \(^ |2/-<*/-< (2iv~\)t h^ |2 -ix* i \r< ppixt')^^^ whcu a quantum jump occurs, only the relative phase and 

the total number of particle changes (see equation ([59]) ). 
X e'^* [N - {d„^ + da^ - 1)] [N - {da, + d^^)] in Figdl and FigHS] we compare the constant loss rate 



{bhha) 






-iXt 



X e'^' [N - {da, + da 



l)][N-{da,+da,)] 



N~2 



X Gi{ai,a2) 



{b'^bh 



CfCle 



-(4i?)+A)t 



-2»xt 



X [TV - {da, + da^ - 1)] [A^ - {da, 
X F2(cri,cr2) 



(136)approximation to the exact numerical result in the case 

of overlapping Rb condensates with large asymmetric two 

2 2ixt\^^'^ body losses considered in Section [Hjn FigfT2l we address 

I ''I ^ / the case of evenly split condensates Na = N^ = N/2 while 

da2 )] in Figlin] we address the case of breathe-together param- 

(137)eters. 

The constant loss rate approximation neglects two ef- 

where we introduced the functions F/3 (cti , cr2 ) and G/3 (cti , cr2 ) ^^^^^: ^he decrease of the loss rate in time as less and less 

particles are in the system, and the change of the ratio 
{Na)/{Ni,) as particles from the b component are lost. In 
the case of FigfT2l where we consider initially {Na) — {Ni,), 
which is the most favorable for squeezing, these two ef- 
fects partially compensates: one tending to degrade and 
the other to improve the squeezing with respect to reality. 
In the case of FigfT51 instead, the two effects sum-up, both 
of them tending to degrade the squeezing with respect to 
reality. Note however that even for such large and com- 
pletely non-symmetric losses, the constant loss rate ap- 
proximation proves to be rather accurate up to the best 



Ff3{ai,a2) 



exp 



Et 



la 



[1- 



-im/3(x+x)tl 



^^ iml3{x + x)[|Ca|2e-^/3(x+x)t + \Cb\^e'0(^-x)t] 

imf3{x-x)t _ Ij 



„m<T2^(™)[g. 



imf3{x - x)[|CaPe-'/5(x+x)* + |Cfc|2e»'3(x-x)t]r 



'lab[l 



_ g-i2/3xtl 



i2/3x[|ape-'/5(x+x)t _,_ |c^|2g^/3(x-x)t] 



(138) 



squeezing time. 
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Approximate analytical solution 
Full numerical simulation 
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Fig. 12. (Top) Spin squeezing with two-body losses in a bi- 
modal Rb condensate as a function of time for symmetrically 
split condensates. Blue solid line: exact numerical simulation 
with 4000 realizations. Red dash-dotted line: analytical solu- 
tion with constant loss rate approximation. (Bottom) Corre- 
sponding total number of particles and (Sz) as a function of 
time. Parameters: Na ^ Nb = W^ , uj = 2n x 42.6 Hz, m=87 
a.m.u., Uaa = 100.44rB, abb = 95.47 rs, Oab = 88.28r_g, rs is 
the Bohr radius, X = 5.367xl0"^s-\ x = 5.412xl0-''s-\ u = 
13.758s~\ Particle losses: K^"^ = 0, A'f * = 119x 10"^^m^s"\ 
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Fig. 13. (Top) Spin squeezing with two-body losses in a bi- 
modal Rb condensate as a function of time in breathe-together 
configuration. Blue solid line: exact numerical simulation with 
4000 realizations. Red dash-dotted line: analytical solution 
with constant loss rate approximation. (Bottom) Correspond- 
ing total number of particles and (Sz) as a function of time. 
Parameters: Na = 7432, Nb = 12568, x = 5.392 x 10"^s"\ 
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